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ABSTRACT 


The  theory  of  binary  coherent  systems  is  generalized  for 
multi-state  components.  The  system  state  ig  defined  to  be  the 
state  of  the  ‘’worst1*  component  in  the  "best"  min  path,  or  „ 
equivalently,  the  state  of  the  ■"best”  component  in  the  *worst" 
min  cut.  All  of  the  results  for  the  binary  case  can  be  com- 
puted for  multi-state  systems  using  the  binary  structure  and 
reliability  function  concepts.  Monotonicity  results  are  now 
valid  with  respect  to  stochastic  ordering  of  component  proba- 
bility vectors.  , 


COHERENT  SYSTEMS  WITH  MULTI-STATE  COMPONENTS 


by 

R.  E.  Barlow 


1.  INTRODUCTION 

The  theory  of  binary  coherent  structures  has  served  as  a unifying 
foundation  for  a mathematical  theory  of  reliability  [1].  Various  generali- 
zations to  multi-state  coherent  structures  have  been  suggested  ( [2 ] , [ 3] , [4 ] ) . 
However,  these  generalizations  have  not  been  fruitful  except  for  very  special 
applications.  We  define  a system  state  function  for  coherent  systems  with 
multi-state  components  and  investigate  its  properties.  Many  results  for  the 
binary  case  have  natural  extensions  in  terms  of  this  system  state  function. 
These  results  also  have  applications  in  fault  tt _e  analysis. 

Suppose  that  we  have  a system  with  components  C = {1,2,  ...,  n}  . Fur- 
thermore, suppose  that  each  component  can  be  in  one  of  m + 1 states, 

{0,1,^,  ...,  m}  , where  0 is  the  failed  state  and  m is  the  maximal  or 

"perfect"  state.  In  addition,  we  have  given  sets  of  components  called  min 

P 

path  sets  ***»  ^ ^ where  C = U Pr  . No  min  path  is  properly 

P r=l  r 

contained  in  any  other  min  patn  set.  The  components,  C , together  with  the 
min  path  sets  define  a coherent  system.  The  system  structure  determines  the 
min  path  sets.  Intuitively,  if  all  components  in  at  least  one  min  path  set 
are  "functioning",  then  the  system  is  "functioning".  This  is  the  set  theo- 
retic definition  of  a coherent  system  [cf.  [5]J.  The  min  path  sets  determine 

a blocking  collection  of  sets  called  the  min  cut  sets  ...»  K^}  . 

k 

Each  min  cut  set  meets  each  min  path  set  and  C = U K . Also,  no  min  cut 

s-1  S 

set  is  properly  contained  in  any  other  min  cut  set.  Intuitively,  if  all  com- 
ponents in  at  least  one  min  cut  set  are  "not  functioning",  then  the  system 
is  "not  functioning".  Either  the  min  path  sets  or  the  min  cut  sets  uniquely 


2 


determine  our  coherent  system. 

Let  z^  = j if  component  i is  in  state  j (0  j m)  so  that 

£ = (z^,Z2,  ...»  z^)  is  the  component  state  vector.  The  specification  and 

determination  of  component  states  will  in  general  depend  on  engineering  and 

system  considerations  which  we  will  not  discuss  here.  Let  capital  be 

m 

a random  variable  and  let  P(Z.  = i)  = p . . >0  where  J"  p.,  = 1 . Since 

1 iJ  - j=0 

p could  be  0 for  some  states,  it  is  not  necessary  that  every  component 
be  capable  of  assuming  every  state.  In  general,  component  states  will  be 
qualitative  measures  as  are  the  concepts  "failed"  and  "functioning".  In 
much  of  what  follows,  it  is  not  necessary  to  confine  component  states  to 
integer  values  {0,1,2,  ...,  m}  . Z^  could,  for  example,  take  values  in 
[0,1]  or  even  negative  values.  However,  we  do  not  pursue  this  generaliza- 
tion here. 

The  performance  level  of  a system,  given  the  component  state  vector 

z^  = ...,  z^)  , will  be  system  dependent  and  it  is  unlikely  that  any 

one  mathematical  definition  of  system  performance  will  be  preferred  above  all 

others.  Hence,  we  concentrate  on  a fundamental,  but  necessarily  limited 

measure  of  system  performance.  If  the  coherent  system  is  a series  system, 

then  we  assign  to  the  system  the  state  of  its  "worst"  component;  i.e.  if  c 

is  the  system  state  function,  then  C,(z)  = Min  z.  . Intuitively,  a series 

l<i_<n 

system  is  no  better  than  its  worst  component.  If  the  coherent  system  is  a 
parallel  system,  then  we  assign  to  the  system  the  state  of  its  "best"  com- 
ponents; i.e.  C( z)  = Max  z . We  will  make  use  of  the  following  well 

l<i£n  1 

known 


1.1  Proposition: 

For  a coherent  system  with  min  path  sets  ...,  Pp)  and  min 

cut  sets  {K^,K2»  ....  K^J  and  any  real  valued  function  f^  , 


r” 


V ’ 


Max  Min  f.  = Min  Max  f. 

l<r<p  ieP  1 l<s<k  ieK 

r s 


[See  [1],  p.  12  and  15.] 


1.1  Definition: 


For  a coherent  system  with  min  path  sets  Pp^  and  min  cut 


sets  {K^,^,  ••••  the  system  state  function  is 


(1.1) 


4(z)  = Max  Min  z = Min  Max  z 

l<r<p  ieP  l<s<k  ieK 

r r s 


Note  that  ^(z)  is  coordinate-wise  nondecreasing.  Intuitively,  t(z)  is 
the  state  of  the  "worst"  component  in  the  "best"  min  path  set  or,  equiva- 
lently, the  state  of  the  "best"  component  in  the  "worst"  min  cut  set. 

With  this  definition  of  system  state,  most  of  the  results  for  binary 
coherent  systems  have  a natural  generalization.  For  example,  suppose  that 
t^  is  the  first  time  that  component  i reaches  state  j starting  in 
state  m , then  the  time  until  the  system  first  reaches  state  j starting 
in  state  m , t , is  easily  seen  to  be 


(1.2) 


t.  *=  Max  Min  t = Min  Max  t 

^ l<r<p  ieP  ^ l<s<k  ieK  ^ 
r s 


The  result  that  redundancy  at  the  component  level  is  better  than  re- 
dundancy at  the  system  level,  [Theorem  2.4,  p.  8 [1]],  also  has  a natural 

generalization.  Lei  z = (z,,z0,  ...,  z ) and  w = (w,,w„,  ...»  w ) be 

iz  n — lz  n 

component  state  vectors  for  components  (1,2,  ...,  n)  and  {1P,2P,  ...,  n^} 
respectively.  Components  i and  i^  may  be  identical  but  operate  inde- 
pendently of  each  other. 


z V w = (z.  V w. , z.  V w0,  ...,  z V w ) 
— — i 1 z z nn 


Define 
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where  z±  V w±  = max  (z^w.,)  . Then 

(1*3)  4(z.  V w)  _>  s(z)  V C(w)  . 

This  is  an  immediate  consequence  of  the  coordinatwise  nondecreasing  property 
of  t . [Definition  1.1.] 
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2.  STOCHASTIC  PROPERTIES  OF  THE  SYSTEM  STATE  FUNCTION 


x.  . = 
13 


J1  1 


f component  i is  in  state  j 


0 otherwise. 


Clearly  £ x.  . = 1 . Let  be  the  usual  coherent  structure  function  as- 
j=0 

sociated  with  min  path  sets  {P,  ,P~,  ....  P } . Then  4>(x,  ,x„ x ) = 1 

1 2 p lm  2m  mm 

if  and  only  if  at  least  one  min  path  set  has  all  of  its  components  in  state  m . 


Also  <t>(x,  .x,,  .....  x ) = 0 if  and  only  if  at  least  one  min  cut  set  has  no 
Im  zm  nm 

components  in  state  m . Thus  $ is  the  coherent  structure  function  which 
recognizes  only  two  states:  the  state  m (or  perfect  state)  and  the  nonper- 

fect set  of  states  less  than  m . 
m 

Let  y±.  = l x±r  and  ^ •••.  Ynj)  • Clearly  - £3+1 


coordlnatewise,  so  that 


Also  <t>(y.j)  = 1 if  a°d  only  if  there  is  at  least  one  min  path  all  of  whose 
components  are  in  state  j or  greater.  Hence  r,(z)  j if  and  only  if 
4>(£j)  = 1 . It  follows  that 


(2.1) 


<K£j)  " ♦ (Xj+1>  = 


if  C(z)  = j 


0 otherwise. 


This  observation  will  be  helpful  for  computing  P(i;(Z)  = j)  since 


(2.2) 


P(C(Z)  = j)  = P[*(Y  ) - <KYj+1)  = 1] 


= E*(Yj)  - E*(Y  x)  . 


In  the  following  we  assume  that  all  components  are  statistically  independent. 
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We  have  proved  the  following 
2.1  Theorem: 

m 

Let  P(X„  = 1)  = p„  and  q„  = £ p^r  . For  a coherent  system  with 

structure  function  <(>  and  reliability  function 

h(£m>  = E<f>(Xlm’X2m’  *•*»  Xnm) 


k 


we  have 


(2.3) 


P(C(Z)  = j)  = h(c^)  - h(£j+1)  0 < j £ m - 1 

P ( C (2)  = rn)  = hCgJ  . 


(2.4)  P(c(Z)  >.  j)  = h(£j)  0^j£m, 

where  £ = (q1J,q2j,  ....  qn;j)  • 

Examp le : 

Consider  the  following  two  terminal  network  representation  for  a three 
component  coherent  system.  Assume  components  can  be  in  any  one  of  three 
states:  0 for  failed,  1 for  marginal  and  2 for  perfect.  The  usual 


FIGURE  2.1.  Two  Terminal  Network, 
reliability  function  is 


h(£2)  * ?12  22  E ^32^ 
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where  p22  u p32  = p22  + p32  - P22P32  * Let  q-ri  = P-n  + P-i?  and  q-o  = P 


il  Kil  Ki2 


*12  Ki2 


for  i < i < 3 . Then 


P (system  is  marginal)  = P(£  = 1) 


= h(£^)  - h(£2) 


= qu(q21  n q31)  - q12(q22  u q32) 


Notation: 


It  will  be  convenient  to  let  p = (p_)  be  the  probability  matrix  cor- 
responding to  component  state  matrix  x = (x„ ) . Let 


(2.5) 


hj(g)  = h(fLj)  “ h(£.j+1)  0 i J 1 ra  “ 1 


and  hm(P)  = 

be  the  probability  that  the  coherent  system  is  in  state  j and  let  h(p)  = 
(hp(p) ,h^(p) , hm(p))  . The  usual  monotonicity  properties  of  binary  co- 

herent systems  have  analogous  monotonicity  properties  with  respect  to  stochastic 

ordering.  Let  a = (ct^.a^,  a^)  be  a probability  vector;  i.e.  0 <_  a < 1 

m 

and  \ a = 1 . 
j=0  J 


Definition: 


a < a iff 


m m 

l ai  1 l a r 

r-j  J r=*j 


for  0 < j < m ; 


i.e.  a is  stochastically  less  than  a . 


Definition: 


E 1 

st 


iff 


(PiO,Pil’  Pim)  - \PiO,Pil’  *•••  Pim) 

st  ' 


for  i _<  i _<  n . Intuitively,  the  better  the  components,  the  better  the 
system  state.  More  precisely,  we  have 


2.2  Proposition: 

* 

If  p < p , then 
st 


(2.6) 


h(p)  £ h|p  j 


Proof: 


To  show 


? ? / *\  / *\ 
l h (p)  <_  l hip)  we  need  only  verify  that  h(£  ) h(£  ) 

r=j  r = r=j  r\=  / J V J / 


■ nr  m i , 

A p p A A 

by  Theorem  2.1  (2.4).  Since  p _<  p implies  q.  = I pir  <_  l pir  = q.. 

= st  = 1J  r=j  ir  r=j  ir  1J 

and  h is  coord inatewise  nondecreasing,  the  result  is  obvious. | 


Generalization  of  the  Moore-Shannon  Theorem 

Moore  and  Shannon  showed  that  binary  coherent  reliability  functions  are 
S-shaped  in  the  sense  that  if  all  components  work  with  probability  p , either 
h(p)  _>  p or  h (p)  <_  p for  all  p , or  there  exists  p^  such  that  h(p)  <_  p 
for  p <_  Pq  while  h(p)  p for  p _>  Pq  • [Cf.  Theorem  5.4  [1].]  This  re- 
sult, comparing  an  arbitrary  binary  coherent  system  reliability  with  a single 
component  reliability,  has  a natural  generalization  with  respect  to  stochastic 


ordering. 


9 


2. 3 Proposition: 


Let  (P-^.P-q.  •••.  Pim)  = (ao.ai»  •••.  «m)  - “ • Assume  h(pQ)  = pQ 

(0  < Pg  < 1)  is  the  fixed  point  for  the  corresponding  binary  coherent  re- 

* 

liability  function.  Let  a = (1  - Pg,0,  . ...  0,Pg)  . Then 


(a)  2—2  implies  h(a)  < a , while 

st  st 

(b)  oi  _>  oi  implies  h(a)  2 2 

st  st 


where  a is  the  probability  matrix  with  identical  rows,  a 


Proof : 


To  show  (a)  we  need  only  verify 


r ? I *\ 

l hr(g)  1 l . 

r-j  r-j 

m j m \ m m ^ 

But  l nr(a)  = h(£.)  = hi  £ ar)  and  £ < l »r  = PQ  for  1 1 j <_ 

r-j  J \r=j  / r=j  r=j 

(m  v ni 

).  a i — I a f°r 

r-j  r/  r-j  r 

0 £ j < m which  in  turn  implies  h(a)  £ . 

st 

(b)  is  proved  similarly. | | 


Proposition  2.3  allows  us  to  compare  arbitrary  coherent  systems  (with 
identical  components)  to  a one  component  system  with  the  same  probability 
vector.  Of  course,  if  no  Pq(0  < Pq  < 1)  such  that  h(Pg)  = Pq  exists, 

then  either  h(a)  £ a or  h(a)  _>  a for  all  ot  . 
st  st 
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y 


A Geometric  Property  of  the  System  State  Distribution 


A fundamental  result  in  coherent  structure  theory  is  that  the  distribu- 
tion of  time  to  first  system  failure  is  IFRA  (for  increasing  failure  rate 
average)  if  component  life  distributions  have  the  same  property.  The  same 
kind  of  result  carries  over  to  the  system  state  distribution. 

Suppose  states  tj,j  + 1,  . . . , m)  are  the  "good"  states  and  the  distri- 


bution of  time  for  a component  to  leave  the  good  states  starting  in  state  m 

,^1/t 


is  IFRA;  i.e.  {PfZ^t)  _>  j 1 ) is  nonincreasing  in  t _>  0 for  fixed  j 
Then 


(PU(t)  1 j]} 


1/t 


is  also  nonincreasing  in  t _>  0 for  fixed  j . This  is  the  so-called  IFRA 
closure  theorem  [Cf.  Theorem  2.6,  [1]].  The  following  is  a corollary  to  the 
IFRA  closure  theorem. 


2.4  Proposition; 

If  (PtZ^Ct)  j is  nonincreasing  in  j _>  0 for  fixed  t , then 

1/3 


{Ptc(t)  > j]} 
is  also  nonincreasing  in  j > 0 for  fixed  t 
Proof : 


PU(t)  >_  j]  = h [ ( t ) ] where 


<?1j  (t)  = P(Z1(t)  > j] 


Define  F (j ) = \ p (t)  and  elsewhere  by  constructing  the  linear  inter- 

r-j  lr 

. — 1 / X 

polate  to  the  points  (j,-log  F^j))  j = 0,1,  m . Then  [F^fx)]  ' 


is  nonincreasing  in  x _>  0 . By  the  IFRA  closure  theorem 


r 


- 


f 


{[hlF1(j),F2(j),  ....  Fn(j)]}1/j 

= [h[q1j(t),q2;j(t) qnj(t)^1/:i 

is  nonincreasing  in  j _>  0 . | | 

Measures  of  Component  Importance 

Various  measures  of  component  importance  have  been  suggested  for  com- 
ponents in  a binary  coherent  system  [6],  [7].  These  have  a natural  general- 
ization to  multi-state  systems.  For  example,  we  say  that  component  i is 
"critical"  to  a coherent  system  at  state  j , if,  with  component  i in  state 
j , the  system  is  in  state  j and  with  component  i not  in  state  j , the 
system  is  not  in  state  j . Let  1^  (p)  be  the  probability  of  this  event. 

We  call  1^  (p)  the  probability  importanae  of  component  i with  respect  to 
system  state  j . 

To  compute  I (p)  , let 

(1ij’^j)  = (yiO’yil yi, j-l,1,yi, j+1*  •••’  yim)  * 

Then 


4,(1  ij’V  - ♦<°ij.Xj+i>  = 1 

if  and  only  if  component  i is  critical  to  system  states  {j,j+l,  ....  m) 
in  the  sense  that  if  component  i is  in  state  j , then  the  system  is  in 
one  of  the  states  {j,j  + 1,  ...,  m}  and  if  component  i is  not  in  state 
j , the  system  is  not  in  one  of  the  states  {j,j  + 1,  ...,  m)  . Hence  com- 


ponent i is  critical  to  system  j if  and  only  if 


L_ 


“ ♦<°ij.r|>]  ■ [<t>(1i,j+l’ij+l^  “ ^i.j+l’^j+i)]  = 1 

for  j = 1,2,  ...,  m - 1 . Hence 


m - 1 . 


Hence 
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V*)  = E[*(1ij’V  - *<W]  - Et4'(1i,j+l’-j+l)  ' <(>(0i,j+1,-j+l>] 


(2.7) 


[hdij.Sj)  - - [h(1i, j+i»aj+i)  - h(0i,j+i*5-j+i)] 


for  j = 1,2,  m - 1 . It  is  easy  to  verify  that 


Iim(I)  = ^^im’V  - *<0im’V] 


I^(P)  = E [ <t> (1  jL1 »— 1>  " • 


il—1 


Note  that  for  the  case  m = 1 , I.  (p)  = I.  (p)  . 

io  = lm  = 

Finally,  the  probability  that  component  i is  in  state  j and  component 
i is  critical  for  system  state  j is  p^I  (p)  . 

First  Passage  Distribution  to  State  j 

We  now  assume  tha"  components  can  "degrade"  through  successively  lower 
states  until,  finally,  total  failure  coincides  with  passage  to  state  0 ; i.e. 


m -*•  m - 1 . 2 -►  1 -*■  0 . 


If  state  transition  times  are  independent  and  exponentially  distributed,  this 
is  called  a "pure  death  process".  In  this  section  we  do  not  allow  transitions 
to  higher  states.  Let  F^  be  the  (continuous)  distribution  of  time  until 
component  i first  reaches  state  j , starting  in  state  m and  F (t)  = 

1 - F (t)  . Let  s(t)  be  the  state  of  the  system  at  time  t , then 

PU(s)  >j;0£s_<t}  = hlF1J(t),F2j(t),  ...,  FnJ  (t)  ] 
by  (2.4)  of  Theorem  2.1. 

We  can  also  compute  the  probability  that  component  i causes  the  system 
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I 

I 


to  enter  state  j , in  the  sense  that  the  system  state  changes  from  state 
j + 1 to  state  j simultaneously  with  a similar  state  transition  by  com- 
ponent i . If  component  first  passage  distributions  to  state  j are  con- 
tinuous, as  we  assume,  then  at  most  one  component  can  be  responsible  for  the 
system  changing  state. 

2.5  Proposition: 

If  first  passage  distributions  of  time  to  go  from  state  m to  state  j 
(0  _<  j _<  m - 1)  , F^  (t)  , are  continuous,  then  the  probability  that  component 
i causes  the  system  to  pass  to  state  j by  time  t is 

t 

J*  {[h(ly,Fj(u))  - h(o±j  ,Fj  (u))  ] 

0 

- [h(lifj+1,Fj+1(u))  - h(o1>J+1,|j+1(u))]}dFlj(u)  . 

where  F (u)  = (F^  (u)  .F^  (u) , ...,  F^(u))  . 

Proof ; 

By  (2.7) 

[h(li  ,F  (u))  - h(oij ,F  (u)) ] 

[h(1i,j+i,Fj+i(u))  - h(°1^ j+1>Ij+1(u)) ] 


is  the  probability  that  component  i is  critical  to  system  state  j + 1 

at  time  u in  the  sense  that  if  component  i is  in  state  j + 1 at  time  ; 

u , the  system  is  in  state  j + 1 at  time  u . Otherwise  the  system  is  not 

in  state  j + 1 at  time  u . If  we  multiply  this  by  the  probability,  j 

dF^j (u)  , that  component  i enters  state  j in  the  interval  (u,u  + du)  , 

__  J 
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and  integrate  we  have  the  required  probability.  Note  that  if  component  i 
is  in  state  j + 1 and  critical  to  system  state  j + 1 at  time  u , then 
it  is  a "worst"  component  in  a "best"  min  path.  Hence  if  component  i then 
leaves  state  j + 1 for  state  j , so  does  the  system  since  component  i 
was  critical. | | 
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